Stabilization of Skyrmion textures by uniaxial distortions in noncentrosymmetric 

cubic helimagnets 
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In cubic noncentrosymmetric ferromagnets uniaxial distortions suppress the helical states and 
stabilize Skyrmion lattices in a broad range of thermodynamical parameters. Using a phenomeno- 
logical theory for modulated and localized states in chiral magnets, the equilibrium parameters of 
the Skyrmion and helical states are derived as functions of the applied magnetic field and induced 
uniaxial anisotropy. 

These results show that due to a combined effect of induced uniaxial anisotropy and an applied 
magnetic field Skyrmion lattices can be formed as thermodynamically stable states in large intervals 
of magnetic field and temperatures in cubic helimagnets, e.g., in intermetallic compounds MnSi, 
FeGe, (Fe,Co)Si. We argue that this mechanism is responsible for the formation of Skyrmion states 
recently observed in thin layers of Feo.5Coo.5Si [X.Z.Yu et al.. Nature 465(2010) 901]. 
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Multidimensional localized and modulated structures 
(Skyrmions) are intensively investigated in many areas 
of physics. [l|, Q In the majority of nonlinear field mod- 
els, Skyrmionic states appear only as dynamic excita- 
tions, but static configurations are generally unstable and 
collapse spontaneously into topological singularities. 
These instabilities can be overcome if the energy func- 
tionals contain (i) contributions with higher-order spa- 
tial derivatives {Skyrme mechanism),^ or (ii) terms lin- 
ear with respect to spatial derivatives (so called Lifshitz 
invariants) Q 



LjdkLi 



(1) 



where L is a vector order parameter (e.g. the magne- 
tization vector M in magnetic materials or the director 
n in chiral liquid crystals), dkLi = dLi/dxk are spatial 
derivatives of the order parameter. 

In condensed matter physics there are no physical in- 
teractions underlying energy contributions with higher- 
order spatial derivatives. [3] On the contrary, the invari- 
ants of type (HJ arise in systems with intrinsic [l^l and 
induced chirality.Q Particularly, in noncentrosymmetric 
magnetic materials such interactions stem from the chi- 
ral part of spin-orbit couplings (Dzyaloshinskii-Moriya 
interactions) . [9] Chiral interactions of typeffl stabilize 
helical |Tlj and Skyrmionic structures |1, |T2| with fixed 
rotation sense (Fig. [1]). Theoretically, isolated Skyrmions 
and Skyrmion lattices have been investigated for several 
classes of noncentrosymmetric systems (e.g. see Refs. 
0, [13 and bibliography in Ref. 0). Contrary to uni- 
axial chiral ferromagnets from Laue classes C„y and 
where thermodynamically stable Skyrmionic states exist 
in a broad range of applied magnetic fields and temper- 
atures, 0, [1^ I13I I Skyrmionic states compete with one- 
dimensionally modulated helices in cubic helimagnets. [1] 
Hence, in the main part of the magnetic phase diagram 



for cubic helimagnets, the ordered state helical structures 
correspond to the global energy minima (such texture 
have been recently observed in (Fe,Co)Si alloys [3] and 
in FeGe fl5'|. Therefore, additional effects are necessary 
to stabilize Skyrmionic states in these svstems.jil. [l3| In 
this paper we demonstrate that uniaxial distortions sup- 
press the helical phases and enable the thermodynamic 
stability of the Skyrmion lattice in a broad range of ap- 
plied magnetic fields. The calculated magnetic phase di- 
agram allows to formulate practical recommendations on 
the possibility to stabilize Skyrmion states at low temper- 
atures in MnSi, FeGe, (Fe,Co)Si and similar intermetallic 
compounds with B20-structure. 

Follow ing the phenomenological theory developed in 
Refs. 0, [Hi we write the magnetic energy density for a 
cubic helimagnet with uniaxial distortions along z-axis 
as 



A (gradM) -M-H + wd 



- km: 
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where A is the exchange stiffness, the second term is the 
Zeeman energy, wd ^ D {A'^^J +A^J'J +A'fJ) =DM-rotM 
is the chiral energy with the Dzyaloshinskii constant D, 
Wa = J2^[B{^iMif + KJAf] includes exchange {B) and 
cubic {Kc) anistropies [11|. The last term in ^ is uni- 
axial anisotropy induced by distortions. 

The Dzyaloshinskii-Moriya energy wd © favours spa- 
tially modulated chiral states where the magnetization 
rotates with a fixed turning sense in the plane perpen- 
dicular to the propagation direction (Fig. [T]). The sign 
and magnitude of the Dzyaloshinskii constant D deter- 
mine the modulation period and the sense of rotation, 
respectively. Thus, in zero magnetic field, H = 0, and 
for zero anisotropics, B — Kc = K — a, flat helix 
forms the magnetic ground state as a single harmonic 
mode with wave number go — D/{2A), where the phase 
angle (f) of the magnetization varies linearly along the 
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FIG. 1: (Color online) Chiral modulated states in noncentrosymmetric cubic magnets: the helicoid, a distorted helix with the 
propagation direction perpendicular to the applied field (a), the conical helix propagates along the applied field (b), an isolated 
Skyrmion (c), a hexagonal Skyrmion lattice (d) with the internal structure of the unit cell (e). 



propagation direction ^, 4>{£) = S,qo.%,^^ Intrinsic cubic 
anisotropy Wa is much weaker than the energy terms in 
Eq. ([2]) and are neglected in further calculations. Its role 
will be discussed to the end of the paper. The solutions 
for one-dimensional modulations include (i) distorted he- 
lices (helicoids) (Fig. [T] (a)) and (ii) conical phases, he- 
lices with the propagation vector along the applied field 
(Fig. [11(b)). 

For the latter state the equilibrium parameters are 
readily derived in analytical form 

cos9^^,^^-^, H^^hJi-^), (3) 

-Ho ^^D \ i^O J 

and the equilibrium energy density Wc — 
-KoM^[H^/{HoHd) - 1]. In Eqs. ^ Ld ^ 2A/D 
is the characteristic length unit of the modulated 
states. In the critical field Ho{K) the conical helix 
flips into the saturated state. The characteristic field 
Hd = D"^ M / {2A) is the flip field for zero distortions. 
The anisotropy value Kq = Hj_/{2M) marks the critical 
value for uniaxial distortions suppressing the conical 
phase in zero field. The analytical solutions for helicoids 
have been derived by Dzyaloshinskii. 

In addition to the helical phase the model Eq. ([2]) 
has solutions for two-dimensional modulated states 
(Skyrmions).[l, To describe Skyrmionic states in 
a magnetic field along the z-axis we introduce spher- 
ical coordinates for the magnetization vector M = 
M (sin 6* cos ■0; sin sin '0; cos 0)) and the cylindrical coor- 
dinates for the spatial variable r — L o{p cos if\ psvaip; z). 
Minimization of energy ([2]) yields rotationally symmetric 
solutions (/3 — ip+TT /2 and 6 = 6{p) derived from equation 



-sm^e- f{e) = 0, 
P 



f{e) = (K/Ko) sin 61 cos 61 + {H/Ha) sin 9, 



(4) 



with boundary conditions, 0{O) = tt, 9{oo) = 0. 

Within a circular cell approximation the equilibrium 
parameters of Skyrmion lattices can be derived by in- 
tegration of Eq. with boundary conditions 9{0) = 
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FIG. 2: (Color online) Equilibrium size Rc of the Skyrmion 
cell and the characteristic radius of the core Ro as functions 
of an applied magnetic field for K/Ko = 0.8. Inset (a) shows 
the evolution of Skyrmion profiles 9{p) with increasing mag- 
netic field, the distributions of the perpendicular magnetiza- 
tion (M^) are sketched in Inset (b): H/Hd = (1, I), 0.2 
(2, II), 0.454 (3, III), 0.6 (4), 0.8 (5) (profiles (1,2) describe 
the magnetization in lattice cells (I, II), while profiles (4,5) - 
isolated Skyrmions. Profile (3) corresponds to the transition 
of the Skyrmion lattice into a system of isolated Skyrmions 
(III). 



TT, 9{R) — and a subsequent minimization of the lat- 
tice energy density Wsl = {2lE?) w{p)pdp with re- 
spect to the cell radius Mathematically similar 
equations arise for Skyrmion states in uniaxial noncen- 
trosymmetric ferromagnets. [l^ In particular, there are 
crystal classes where the Dzyaloshinskii-Moriya energy 
is described by Lifshitz invariants with gradients only 
along the directions perpendicular to the unaxial axis. 



e.g. for D2d classes, wd = D{K 



(x) 
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This restriction of modulations to two dimensions proved 
to be crucial for the thermodynamical stability of the 
Skyrmion states. [13, [13: [iBl In the high symmetry cu- 
bic helimagnets the chiral energy wd = -DM • rotM ^ 
energetically favours the helical phases compared to the 
Skyrmion states. Q 

From the numerical investigation of Eq. we now 



3 



show that a sufficiently strong magnetic anisotropy K 
stabihzes Skyrmionic textures in apphed magnetic fields. 
Typical solutions for 9{p) are plotted in the inset of Fig. 
[21 At zero field the Skyrmion cells have a smooth dis- 
tribution of the magnetization (profile (1) in Fig. [21 
Inset). An increasing magnetic field gradually squeezes 
the Skyrmion core (profiles (2), (3)) and transforms the 
lattice into a system of isolated Skyrmions (profiles (4) , 
(5)). The Skyrmion core size can be introduced in a 
manner commonly used for magnetic domain walls, Q 
^0 = Po + &o/{dO/dp)o where (po:^o) is the infiection 
point of profile 9{p) (Fig. [2]). The plots in Fig. [2l demon- 
strate a progressive localization of the Skyrmion core Rq 
accompanied by the expansion of the lattice cell size Rc 
with increasing H. 

The functional ([2]) includes two independent control 
parameters, K/Kq and H/Hd- The calculated mag- 
netic phase diagram in these variables (Fig. [3]) provides 
a comprehensive analysis of model ([2]). For K = Q 
the conical phase is the globally stable state from zero 
field to the saturation field {Q < H < while 
metastable solutions for the Skyrmion lattices and he- 
licoids exist below the critical fields Hc/Hd = 0.8132 
and Hb/Hd = = 0.6168, correspondingly. d 

The conical phase is the global minimum of the system 
within the region {a — A — B — d) and transforms discon- 
tinuously into the Skyrmion {{A — S)-line) and helicoid 
phases {{a — A)-lme) . A sufficiently strong K suppresses 
the conical states, and only modulations with the propa- 
gation vectors perpendicular to the applied field can exist 
(helicoids and Skyrmion lattices). The Skyrmion states 
are thermodynamically stable within a curvilinear trian- 
gle (A-B-D) with vertices (A) = (0.050,0.2158), (B) = 
(0.3628,0.6374), and (D) = (1.9004,0.10)) (Fig. [3]). The 
solutions for helicoids exist within area (a — b — D — e). 
For a certain value of magnetic field the helicoid is trans- 
formed into the homogeneous state at the critical line 
(b - D - e) D = (7r/4)VA^[^/TTl^ + arcsinh(l/VI^)] 
where v = H/{2KM).[12] Point (e) designates the criti- 
cal value of the anisotropy K for the suppression of the 
modulated states in zero field {Ke/Ko — 7r^/4 — 2.467). 

The critical points A, B, D separate the phase dia- 
gram (Fig. [31) into four distinct regions. (I) In the low 
anisotropy region {K < Ka = 0.05^0) only helical states 
are realized as thermodynamically stable phases. (II) For 
Ka < K < Kb = 0.363i4ro the Skyrmion lattice becomes 
absolutely stable in a certain range of the applied field. 
The Skyrmionic phase is separated from the helicoidal 
and conical states by first-order transitions. (Inset(c)). 
(Ill) For Kb < K < Kd = 1.90^0 the magnetization 
curve includes a first-order transition between the heli- 
coid and the Skyrmion lattice and the second-order tran- 
sition of the Skyrmion phase into the saturated state (In- 
set (d)). (IV) Finally for { Kd < K < K^) the helicoids 
are thermodynamically stable in the whole region where 



H/K 




FIG. 3: (Color online). Magnetic phase diagram of the so- 
lutions for model ([2]|. Filled areas designate the regions of 
global stability [A — B — D) and meta-stablity (a — 6 — c — e) 
of Skyrmion lattices. These regions transform by first-order 
processes into the helocoids [A — D) or conical phases {A — B). 
The solutions for helicoids exist within the area (a — b — D — e). 
Inset (a) gives the differences AE between the energies for 
the Skyrmion lattice and the conical phase (solid) and the 
distorted and conical helices (dashed) as functions of the ap- 
plied field. Insets (b), (c), (d) show magnetization curves 
for different values of uniaxial distortions: K/Kq = 0.04 (b), 
0.4 (c), and 1.6 (d). Solid lines indicate the globally stable 
phases, metastable states are shown by green thin (conical), 
red dashed (helicoids), blue dotted (Skyrmions) lines. 



the modulated states exist. 

The equilibrium energies of the conical {Wc), heli- 
coidal (Wh), and Skyrmion {Wsl) phases plotted as 
functions AEsi^(^h) = Wsl(h) ~ versus the applied 
field (Inset (a) of Fig. ^ help to elucidate the physical 
mechanisms that lead to the formation of the different 
modulated states. The conical phase has a simple struc- 
ture with its single-harmonic rotation of the magnetiza- 
tion component M± — Msm9 Q. For zero K, i.e., for 
an isotropic helimagnct this provides a larger reduction 
of the energy density in an external field {wd oc —DMf) 
than is possible for the (anharmonic) modulations in the 
alternative phases of helicoids and Skyrmion lattices. An 
increasing uniaxial anisotropy K > gradually decreases 
Q and the chiral energy contribution in the energy 
of the conical phase. Correspondingly, the conical phase 
becomes unstable with respect to the competing modu- 
lated states (Fig. [31). 

The energetic advantage of Skyrmion states is due 
to rotation of the magnetization in two dimensions. 
This (double-twist) grants a larger reduction of the 
Dzyaloshinskii-Moriya energy than a single-direction ro- 
tation in helical phases. Thus, the double-twist yields a 
lower energy density in the Skyrmion cores compared to 
helical states. On the other hand, the incompatibility of 
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spin configurations near the edges of the hexagonal cells 
leads to an excess of the energy density in this region. Q 
The analysis shows that at zero field this energy cost 
outweighs the energy gain in the Skyrmion core. An in- 
creasing external magnetic field anti-parallel to the mag- 
netization in the Skyrmion center gradually decreases the 
total energy by suppressing the energy cost near the wall- 
like structure surrounding the Skyrmion cores with the 
shape of a honeycomb (Fig. [21 Inset). At a finite field as 
marked by the A — D line in Fig. [3l the Skyrmion lat- 
tice has lower energy than the alternative helical states. 
Because the topology of the helix and the Skyrmion lat- 
tices are different, this field-driven transition has to take 
place by a first-order process. Thus, in uniaxially dis- 
torted cubic helimagnets the thermodynamical stability of 
Skyrmion lattices is reached as a combined effect of ap- 
plied magnetic field, that causes the localization of the cell 
core, and the uniaxial anisotropy to suppress the alterna- 
tive conical states. 

Skyrmionic states now have been observed in nanolay- 
ers of Feo.5Coo.5Si To explain the stability of the 

Skyrmion phase in this system, we argue that surface- 
induced uniaxial anisotropy in this system suppresses the 
cone phase and stabilizes skyrmion and helicoid mod- 
ulations in crystal plates that are thin enough. The 
magnetic transformation under field reported for the 
Feo.5Coo.5Si films show a first-order process from heli- 
coids at low field into the dense Skyrmion phase. At 
high fields isolated Skyrmions are set free and form dis- 
ordered ensembles. This sequence of magnetization pro- 
cesses corresponds to the calculated behavior for systems 
with intermediate uniaxial anisotropy. Fig. 3(d). As seen 
in Ref. [2^, there is a type of melting of the ordered 
Skyrmion lattice at higher fields and temperatures where 
the free Skyrmions form weakly coupled disordered ag- 
gregates. This observations corresponds to the existence 
region of free Skyrmions in the magnetic phase diagram 
above the stability range for dense stable Skyrmion lat- 
tices, which is a feature of the generic phase diagrams 
for Skyrmion phases calculated in Ref. [13[. Thus, the 
experimental observations reported in Ref. 22], are in 
close agreement with theoretical predictions on the be- 
havior of Skyrmionic phases, as composed of particle-like 
radial objects (in two spatial dimensions), that remain 
intrinsically stable beyond the existence range of lattice- 
like condensed phases. 0, [13] 

For MnSi, earlier experiments and analysis [l8| of 
magnetoelastic couplings allow a quantitative estimate 
showing that the effects predicted here can be achieved 
in experiment. The magnetoelastic coupling with uni- 
axial strains Uzz is given by Wme — buzz (Mz/Ms)'^ , 
where Mg — 50.9 A/m is the saturation magnetization 
[igf and b = 7.4 GPa is a magnetoelastic coefficient de- 
rived from the magnetostriction data in Ref. [iTl Using 
exchange constant A = 0.11 pJ/m, as estimated from 
the spin- wave stiffness reported in Ref. [2^, and D — 



2qoA = 0.86 /xJ/m^ for MnSi we have Kq ~ 17 kJ/m^ 
and a dimensionless scale b/Ko ~ 44 for the induced 
anisotropy. Thus, a modest strain Uzz = 0.0024 is suf- 
ficient to reach an induced anisotropy K/Kq = 0.1 well 
within the region for stable Skyrmion lattices in the phase 
diagram Fig. [3] This strain corresponds to a tensile stress 
azz— 680 MPa for MnSi by using the elastic constant 
cii— 283 GPa.^] The rather low uniaxial stress neces- 
sary to stabilize the Skyrmion lattice is particularly rele- 
vant for pressure experiments with a uniaxial disbalance 
of the applied stresses, but it could also be achieved in 
epitaxial films. 

Finally, intrinsic anisotropy Wa Q also makes a con- 
tribution to stabilize Skyrmionic states as both exchange 
(B) and cubic (Kc) anisotropics violate the ideal spin 
configuration of the cone ([3]) and increases the en- 
ergy of this phase [13]. This favours Skyrmion lat- 
tices for certain directions of the applied field. Partic- 
ularly, in MnSi {B < 0) for Skyrmion lattices oriented 
along [001] type axes AEsl{H/ H^) = at critical point 
(Bcr = OAOA/D^,Hcr = OAOSHd). Thus, for \B\ > B^r 
Skyrmion lattices are globally stable in a certain interval 
of magnetic fields around Her ■ 

In conclusion, we have shown that in cubic helimagnets 
uniaxial distortions effectively suppress helical states and 
stabilize the Skyrmion states in a broad range of the ap- 
plied fields. 
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